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This paper proposes a new radial basis function neural network maximum 
power point tracking controller based on a differential evolution algorithm for 
machine side converter of permanent magnet synchronous generator wind 
turbine under variable wind speed. Direct axis stator current control methods of 
permanent magnet synchronous machine are reviewed shortly. A combined 
radial basis function neural network-based network maximum power point 
tracking method and d axis stator current control techniques including zero d 
axis stator current, unity power factor, and constant stator flux-linkage have 
been implemented to control the machine side converter of permanent magnet 
synchronous generator wind turbine. The dynamic performance of the 
proposed approach is assessed under different operating conditions through a 
simulation model based on MATLAB. It has been seen that the radial basis 
function neural network controller can not only track well the maximum power 


generator point but also can be reduced costly. 
Radial basis function neural 
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1. INTRODUCTION 

Nowadays, wind energy is popular clean energy in the world due to its own merit are pollution 
freezes, cost reduction and high reliability [1], [2]. Wind energy is interested in both industrial and academic 
fields [3] because of technological advancements and support programs of the government in order to obtain 
green energy strategies for sustainable development [4] 

Generally, the main wind-energy conversion system consists of a wind alternator, a gearbox, a 
transformer and a back-to-back (BTB) converter [5], [6]. Compared to other types, the permanent magnet 
synchronous generator (PMSG) is popular in current uses because of high performance, high reliability, and 
gearless configuration [7] shown in Figure 1. Accordingly, objective of improving performance, efficiency 
and reduce cost, the variable speed wind turbine is needed to operate at the MPPT under variable wind speed 
[8]. The drawback of the conventional MPPT method such as tip speed ratio (TSR) [9] are errors in 
processing measurement, cost and difficulty in implementation [10]. 
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Figure 1. Configuration of PMSG wind turbine 


Likewise, the demerit of optimal controller [11], [12] are the saturation effect and the disturbance on 
the input voltage of the converter. Besides, sliding mode control [13], nonlinear control [14] and fuzzy 
control [15] or robustness control [16] are popularly used due to its highly efficient but its disadvantage is 
failed to track maximum power point (MPP) in the rapidly changing conditions. In the same way, the 
drawback of the maximum electrical power tracking method and a maximum mechanical power tracking 
method [17] or the extreme seeking method [18] is depended on the parameters of the system and maybe 
entrapped in the global maximum. Among the adaptive hybrid intelligent control [19] and the expert control 
system [20], the DE algorithm [21] has received the most attention because of the best effective [22] and 
exact tracking maximum power point [23]. Until now, the DE algorithm has never been implemented in 
PMSG wind turbines. Further, the Radial Basis Function Networks (RBFNNs) [24], [25] based MPPT 
method is proposed to design [26] and implement [27] for PMSG wind turbine because of high performance, 
accuracy, and convergence. Moreover, the d-axis stator current control techniques [28] such as zero d-axis 
stator current (ZDC), unity power factor (UPF) and constant stator flux-linkage (CSFL) have been interested 
in improving performance and reduced cost because of hardware less. Currently, these d-axis stator current 
control techniques have only been applied to control for the motor. In this paper, the study applies these 
techniques with the correction to control for MSC of PMSG. Based on the achieved results, the main 
contribution of this work can be outlined as follows; i) DE algorithm based MPPT is proposed for PMSG 
wind turbine; ii) design and implementation of RBFNN for the zero-approximation error of the DE algorithm 
based MPPT; iii) combined radial basis function neural network based MPPT method and d axis stator 
current control techniques (ZDC, UPF and CSFL) have been successfully applied for the control of the MSC 
of PMSG wind turbine. 

This article is organized as follows: section 2 interprets in detail of wind turbine conversion power. 
Section 3 describes the MPPT method. Section 4 shows the stator current control techniques. Section 5 presents 
and explains the simulation results. Finally, conclusions. 


2. DYNAMIC MODELING OF PMSG WIND TURBINE SYSTEM 
2.1. Wind turbine model 
The mechanical output power of the wind turbine is given by (1), 


p= ~ orR?C, (a, PJV (1) 


~ 2 


where R is the turbine radius (m), p is the air density (kg/m°), vw is the wind speed (m/sec), Cp is the power 
coefficient function, a is the tip speed ratio and £ is the blade pitch angle. 
The tip speed ratio and and power coefficient are expressed, 


a =R (2) 


where c; represent angular speed of rotor 


c5 


c. -D 
C (a, p) = aT C3ß — c4)e % + cea (3) 


where the approximated coefficient values cı — c6 are given in Table 3. 


gi ( 1 2: 4) 


a+0.088 1+8? 


The relationship between turbine power and wind speed under variable wind speed demonstrates in 
Figure 2. The turbine power at a certain wind speed is possible at a maximum, which is called optimum wind 
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speed. The blue curve is described as the optimal tracking curve. The turbine is required to operate at an optimal 
tip speed ratio for maximum power. It can be done by driving the turbine's rotational speed in order that it 
continuously rotates at the optimum speed [2]. The operational zone within a wind speed range that is restricted 
between connected wind speed (vw cut-in) and disconnected wind speed (vw cut-out) affect to the captured wind 
power. Otherwise, wind turbines are required to stop operating above connected wind speed (vw cut-in) or below 
disconnected wind speed (1, cut-out) because of protection conditions. The turbine has to stop for protection of 
the wind turbine as well as generator if out of this range. The wind turbine's rated power (Prata) is determined 
by the wind speed (Varea). Then, it can be divided into 4 main regions and classified as follows [9]; i) wind 
turbines must be stopped and disconnected from the grid to avoid the generator implements in region 1 and 4, 
which is below (vw cut-in) and above (vw cut-out); ii) the second region is in between (Vv, cut-in) and (vw rated) 
which a wind turbines controller implements the MPPT method below rated wind speed to achieve the optimal 
power during the variable wind speed; iii) the third region is in between (vw rated) and (vw cut-out) where the 
pitch controller is used to limit the mechanical power generation and to keep wind turbines in safe operation. 


MPPT region 


Maximum power 


Turbine power (KW) 
Parking region 
Parking region 


cut-in 


Wind speed (m/s) 


Figure 2. Wind energy conversion system operating regions 


2.2. Permanent magnet synchronous generator model 
The equivalent circuit model of the PMSG is depicted in Figure 3. In the dq reference frame, the 
stator voltage of PMSG wind turbine is given by [3]. In this article presents surface mounted alternator. 


_ i dids y 
Uds = Rslas + Las dt T wr Lqsiqs 


digs (5) 
ia Lastds@r + Arwr 


Uas = R.i 
qs s dt 


qs + Las 


Where uas and ug; are the stator terminal voltages on d, q axis, respectively, ias and igs are stator 
current on d, q axis, Las and Lgs are the synchronous inductances of the generator on d, q axis, Rs is stator 
resistance, œ, is the electrical angular velocity of the rotor, 4, is the amplitude of the flux linkage. The 
generator electromagnetic torque, the active and reactive powers at generator terminals are given 


To = Ê piggy (6) 


lol = $v ha + Zug eal (7) 


Where T; is the torque, p is the number of pole pairs 


Lis R: Lpi Lis ip R OLadis Le 


(a) (b) 


Figure 3. Equivalent circuit model of the PMSG: (a) d-axis circuit; (b) q-axis circuit 
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3. MPPT 

In order to increase energy conversion efficiency in varying the speed, WECS need peak power 
extraction because it is important. The MPPT control tries to achieve the maximum possible power at a vw 
certain value of the wind. The maximum power point (MPP) trajectory depends on the wind speed. The 
operational zone of the MPPT control is a period from cut-in wind speed to rated wind speed, as seen in the 
graph of this curve in Figure 2. 


3.1. Conventional MPPT method 

In past years, many MPPT methods have been studied and improved [2], [3]. There are many 
different merits and demerits in each of them. In a way that the most popular method is employed the tip 
speed ratio (TSR) technique [17] that uses (2). Thus, this control technique's own flaw is that it necessitates 
the use of a wind-speed measurement gadget and the pre-determined value of the best tip speed ratio for 
converting wind velocity measurements into their corresponding optimal speed reference [18]. This also 
increases the cost of the system [19]. 


3.2. Differential evolution-based MPPT 

Since 1995, DE was arguably one of the first technical present by R. Storn and K. V. Price for 
global optimization over continuous search space [21], in which its configuration was optimal functions in a 
continuous N-dimensional region. Each individual in the population is an N-dimensional vector that 
represents a problem solution. DE depends on taking the distinction vector between two types and including 
a scaled sort of the distinction vector to a third individual to make another applicant arrangement [22]. 

The process creates a new candidate solution such as a mutant vector is created by joining three 
arbitrarily chose vectors from the number of inhabitants in vectors barring the objective vector. This 
consolidating cycle of three haphazardly chose vectors to shape the mutant vector and a multiplier which is 
the fundamental boundary of the DE calculation [23]. This results in a mutant that might be accepted into the 
population as a new candidate solution. There are two generally utilized hybrid techniques in DE: binomial 
hybrid and exponential hybrid. The fittest vector in each pair is kept for the next DE generation, and the least 
fit is discarded. The basic DE algorithm for an N-dimensional issue is depicted in Table 1 and these 
parameters are simulated in MATLAB. The flowchart for DE algorithm is illustrated in Figure 4. 


Start 


Initialization of DE’s parameters 


v 
Input Vw = rand([3,12],120) 


Y 
Initialization of DE’s initial population 
Wri = Wrmin + rand(Wrmax z Wrmin) 


Y 
Check for limit of variable (wr) 
SE y —___ Table 1. The parameters of DE 

Initialization of an initial equilibrium state Description Value 

Wel, Wel, Wis,..3 FXLEX2FX3,.. Particle number of a generation, Npop 30 

- Maximum number of generations, itermax 50 

Evaluate the target function and update the Crossover probability 0.2 

new equilibrium state Scaling factor lower bound 0.2 

Yy Scaling factor upper bound 0.8 

Update and calculate new variable (wr-new) 


check stop condition 
Pop = Wrmin + rand (Wemax — Wrmin 


Stop 


Figure 4. Flowchart for differential evolution 


3.3. Proposed the RBFNN based MPPT method 

In this section is proposed a model of straightforward neural network architecture, which has only 
one hidden layer and is utilized in many application areas. It is named Radial Basis Function Neural 
Networks (RBFNN). An RBFNN based on MPPT method is designed to enable controlling MPPT without 
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knowing the air density or the mechanical parameters of the turbine. The condition is justified because the 
goal of this study is to track the maximum power point and the reference model provided by (1) justifies the 
demand [24], [25]. The network is trained offline using data acquisition which is implemented by using DE 
that as illustrated in Figure 5. 


Maximum Power 
© 


Pm (p.u.) 


0 0.5 1 1.5 2 2.5 3 3.5 4 
wr (rad/s) 


Figure 5. Proposed RBFNN training based maximum power curve of DE algorithm 


The interpolation issues in a multidimensional space depending on many centers as data points were 
solved by RBFNN, which was first present by Powell. Its significant characteristic is fast and in a network, a 
linear learning algorithm capable of showing complex non-linear charts are used. Moreover, the 
generalization characteristics of RBFNN are also to improve [26]. This neural network designs as following 
steps, firstly we run the initial program with the DE algorithm based MPPT and make a detailed look up table 
of input and output data. Secondly, implement to train the RBFNN based this data and using this neural 
network to replace old MPPT block. The RBFNN has three layers. The first layer is input player that has only 
corresponded the input signal. It's possible that there are multiple variables, each of which corresponds to a 
separate neuron. The input layer neurons’ outputs are used to deliver values to each of the hidden layer 
neurons’ inputs. The structure of the RBFNN shows in Figure 6 (a) to realize data processing duties, which 
contains an input, an output, and a hidden layer. The (8) represents the net input and output of the input layer 
and Figure 6 (b) shows the block diagram of the MPPT algorithm in Matlab. 


net; = x} } (8) 
yi = f (net) = (x7) j=1..n;i=1,2 


Where x’; is the input layer, the net sum of nodes is represented by f(x’; ) which is fed to the hidden layer. 

The second layer is a hidden layer that is nonlinear and uses Gaussian functions. It can contain a 
large number of neurons. This layer's selection of neurons is a difficult undertaking. For the selection of 
neurons in the hidden layer, we proposed a data structure-preserving criterion technique. Each neuron in the 
hidden layer is made up of an RBF centered at a point, depending on the complexity of the input/output 
variables. The basic functions, which are dependent on the parameters {w’);, b/;} and input activations {x/;} 
determine the hidden unit activations. For each dimension, the RBF function's spread (radius) may be 
different. Figures 6 (c) and 6 (d) shows the structure of Layer 1. The best structure for an RBFNN is 
determined by the application area. The network's training process determines the network's centers and 
spreads. The hidden layer's net input and output are shown. 


2 1,,1 
net; = Lwii¥1) 


yi, = f? (wit; + bj) =o 


(9) 
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Figure 6. Design and implement RBFNN controller, (a) block diagram of MPPT method using RBFNN,; 
(b) RBFNN block in MATLAB; (c) the structure of RBFNN in MATLAB; (d) the structure of RBFNN in 
MATLAB 


Where w77; is the weights which connect the input layer and hidden layer, b’; is the bias terms of the 
hidden layer. Thirdly, The Gaussian outputs are combined by the linear layer. During training, only the 
weights between the hidden and output layers are modified and calculated using Moore-Penrose generalized 
pseudo-inverse [27]. 


net? = why? 
j (10) 


yi = f (wiyh +b) 


Where w?;j is the weights which connect the hidden layer and output layer, b?;j is the bias terms of 
the output layer. Today RBFNN research has been continuously increasing and becoming more ubiquitous in 
several application areas. RBFNN broadly consists of three layers feed-forward neural network [24], [25]. 
The three-layer RBFNN employed in this section has been demonstrated to be a universal approximation. 
[26]. The RBFNN's overall input-output mapping is given by, 


<r = b+ yy exp (-Ear) (11) 


j 
where x = [vy] are the input vector, C; € Rn is the center of the jm RBF units in the hidden layer, n is the 


number of RBF units, b and vj are the bias term and the weight between hidden and y = [w,] output layers 
respectively of the RBFNN that represents the wind turbine rotational speed. 
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4. d-AXIS STATOR CURRENT CONTROL TECHNIQUES 
4.1. ZDC 

The ZDC is the most widely utilized control technique applied to industry because of its simplicity 
[2]. In order to implement this control technique, d-axis stator current is set to zero. So, the relationship 
between torque and current can be linearized. 


4.2. UPF 

Under this control law, the power factor angle is always kept at unity. This reason would lead to a 
cost effective solution for the BTB converter. This is one of the most important issues for a megawatt-scale 
wind energy conversion system [3]. (12) shown this implementation of this control law by setting the reactive 
power equal zero and illustrated in Figure 7, as shown in (5) substitutes into (12). 


Uasltgs m Ugslas = 0 (12) 


q-axis 


ZDCYV i. 


‘qsZDC 


d-axis stator current control 


igs = — 22+ |(2) = (igs)” (13) 


2Ls 


In (13) shows that the current of g-axis stator should satisfy the following constraint 


l SST (14) 


4.3. CSFL 

In order to overcome increasing stator flux linkage (saturation of the stator yoke) when the torque 
value is increased, the CSFL must be applied. One of the most characteristics of the CSFL control is high 
steady-state efficiency, good capacity factor, and tiny required power converter capacity. In a similar way 
with implementing other d-axis stator current control techniques, the amplitude of the stator flux vector is 
kept constant and equal to the permanent magnet flux shown as (15). This control technique will be done and 
it is illustrated in Figure 7 [28]. 


n : 2 
lAs| = J (Lsias + Ar)? t (Lsiqs) = Ar (15) 


: Ar AN? >42 
ias = -2+ JE) - i) (16) 
In (16) is existed, the current of g-axis stator should satisfy the following constraint 


< (17) 


igs ST 


In this case, the analysis can be summarized as the main feature of ZDC control is high cost under wind speed 
variations due to poor efficiency and the high reactive power: 
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— The main feature of UPF control reduced the cost of the power circuit because of its small size, though it 
is high performance under wind speed variations. This is one of the significant considerations for 
megawatt-level wind turbine capacity; 

— A proposed technique for controlling the MSC of vector-controlled PMSG wind turbine is CSFL control 
because of good characteristic and high power factor. 


5. SIMULATION RESULTS 

In this paper, The RBFNN controller based MPPT method applied to PMSG has been implemented 
in MATLAB under variable wind speed. In order to evaluate and analyze the proposed solution's 
performance, RBFNN controller based on the average and standard deviation (SD) in Table 2 that is 
relationship between the turbine powers at variable rotor speed and the net sum of nodes is 10. When the 
error value reaches the permitted value, the training ends (performance = 0.00001) or lessons reach the 
allowable value (epoch=100). The process of network training is depicted in Figure 8. The network training 
chart presents network-training process, which occurs 100 times before stopping. Simultaneously, the finest 
training results will be 0.00015265 using the RBFNN block takes the place of the old block of the MPPT 
algorithm and is depicted in Figure 7. 


Table 2. Data of RBFNN 


Signal 1 2 3 4 5 6 7 8 9 10 
Average 2.08e- 1.0le+1 2.43e+1 4.07e+1 1.35e- 8.06e- 2.41e- 1.4e-19  8.05e-1 3.72e- 
14 3 14 11 1 
RE SD 0 1.3e+1  3.85e+1  1.62e+2 Oe+0 5.65e- 1.03e- 3.64e- 1.79e-3 Oe+0 
14 22 36 
Signal 11 12 13 14 15 16 17 18 19 20 
Average 3.45e- 3.1e-1 7.12e+1 2.33e+0 5.3e-2 1.38e-6  3.7e-11 9.07e- 3.33e+0 0 
10 20 
RBFNN $D 495e- 4.52e-1 5.39e-1  1.99e+0  6.Je- 1.21e7 178e- 6.8e-23 1.07e1 0 
12 12 12 


Best Training Performance is 0.00015265 at epoch 100 


Mean Squared Error (mse) 


10€ 1 1 1 1 L L 1 1 1 
o 10 20 30 40 50 60 70 80 90 


100 Epochs 


Figure 8. Network training graph 


The simulation results in Matlab, which assess and compare the strong points and weak points of the 
proposed stator current control techniques and RBFNN based on the MPPT method. The effectiveness of 
stator current control techniques is considered under wind speed variation. Tables 3 and 4 list the PMSG 
wind turbine's parameters. Figure 9 (a) is depicted the variable wind speed from 3 m/ses to 12 m/sec over a 
period of 50 sec. In order to produce maximum power and be avoided demagnetization of the permanent 
magnet, the d axis stator current is set up at zero. The simulation results of the stator current on d axis (iaszpc, 
idsUPF, İascsFL) in three control techniques are illustrated in Figure 9 (b). In Figure 9 (b) shows that the highest 
amplitude value of stator current on d axis is the UPF control technique, oppositely is ZDC control technique. 
The disadvantage of the UPF control technique must select wire with large cross-sectional for MSC because 
of the high cost. On the other hand, the stator current on q axis (igszpc iqsuPF, iqscsrFL) Of three control 
techniques are same value (igszpc = iqsupr = iqscsrL) under wind speed variation and illustrated in Figure 9 (c). 
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The stator fluxes (Azpc, Auer, Acsrr) are illustrated in Figure 9 (d). The stator flux of the CSFL control 
technique always is constant while the stator flux of the ZDC control technique is the highest worth and the 
stator flux of the UPF control technique is the lowest worth under different wind speeds. 


Table 3. Turbine parameters [2] 


Description Parameter Value 
Transparent power Sn 2.2MW 
Current In 2606A 
Voltage Un 690V 
DC voltage Vue 1200 V 
Rotating speed Om 2.355 rad/s Table 4. Generator parameters [2] 
Number of poles Zp 26 N.m PPE 
Depiction Parameter Value 
Moment Te 934.1 kNm 
A . Rated power Py 2MW 
Viscous damping B 0.0041 Nms : 
: Rotating speed On 2.355 rad/s 
Stator phase resistance R; 0.8e-3 Q - i 2 
s Alternator rotor inertia Tin 0.25e3 kg.m 
Stator phase inductance L; 1.67 mH 
; Blades length R 37.1 m 

Flux linkage Ar 9.18 Wb wind speed v, 12m/sec 
Inertia of turbine rotor Ja 0.5e6 kg.m? p £ 

ci 0.5176 

C2 116 
The approximate C3 0.4 
coefficient values C4 5 

C5 21 

C6 0.0068 


1 aszpc| 
1 ásUPF 


"ascsril 


Windspeed (m/s) 


5 E 

d-axis stator current (KA) 
+ b 
in ety 


0 5 10 15 20 25 30 35 40 45 50 “0 5 10 15 20 25 30 35 40 45 50 
Time (s) Time (s) 


(a) (b) 


rent (KA) 


stator flux (Wb) 


Time (s) Time (s) 


(c) (d) 


Figure 9. Wind energy conversion system: (a) wind speed; (b) d-axis stator current; (c) g-axis stator current; 
(d) stator flux 


Figure 10 shows the active power (Pzpc, Purr, PcsrL) and reactive power (Qzpc, Quer QcsrL) of 
three control methods, in that the active power is approximately the same value (Pzpc = Pupr = Pcsrt) 
because of according to (6) torque or active power are proportional with g-axis stator current (as mentioned 
above, the Simulink results in Figure 9 (c) shows that the g-axis stator current in three control methods is of 
the same value). Furthermore, the reactive power in the UPF control method is always zero, the highest value 
is reactive power in the ZDC control method. In Figure 11. The highest value of apparent power is the ZDC 
control technique at 12 m/sec and this is its drawback due to a high cost. The opposite is the UPF control 
technique, which has the lowest value. The apparent power of the ZDC control technique continues the 
highest while the UPF control technique and CSFL control technique are the same value over a period of 8 
m/sec to 11 m/sec. Finally, the apparent power of the three techniques is the same value at a wind speed 
below 8 m/sec. The opposite of the highest value of the capacity factor in the UPF technique (@ypr) is the 
ZDC control technique (8zpc) while the power factor of the CSFL control technique (@csrz) has the middle 
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value. It is depicted in Figure 12. According to Figure 13, the stator current of three control techniques (iszpc, 
isupr, iscsrL) has two instances. Over a period of 8 m/sec to 11 m/sec, the highest value of stator current is the 
UPF control technique (isupr) and this is its drawback due to a high cost (it should choose the wire with a 
large cross-section for MSC). The opposite is the ZDC control technique (i;zpc), which has the lowest value. 
When the wind speed is less than 8 m/sec, the stator current of three control techniques (iszpc supr, iscsFL) 
has the same value. 
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Figure 10. Active power and reactive power Figure 11. Apparent power 
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Figure 12. Power factor Figure 13. RMS stator current 


6. CONCLUSION 

In this paper, a new radial basis function neural network maximum power point tracking controller 
based on differential evolution algorithm for machine side converter of permanent magnet synchronous 
alternator is presented. The proposed controller was designed and applied successfully to permanent magnet 
synchronous generator wind turbines. The simulation results show that the controller has high accuracy and 
stability in tracking the maximum power point and increases the power of permanent magnet synchronous 
generator wind turbines. In addition, three d axis stator current, zero d axis stator current, unity power factor, 
and constant stator flux-linkage have been studied. Generally, the constant stator flux-linkage method can be 
reduced costly without the hardware. In future work, a low-cost processor-in-the-loop platform will be 
applied to experimental verification of the proposed approach. 
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